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Abstract 


Commercial RF and microwave CAD packages have evolved to the point 
where predictive simulation is possible. When properly applied, CAD can 
reduce design cycle-time and provide insight that bench techniques alone 
cannot. Virtual prototyping and statistical optimization are possible, 
precluding the necessity of time-consuming pilot runs. Unfortunately, many 
CAD users are unfamiliar with rigorous modeling and characterization 
techniques. The result is often user frustration and alienation. Ultimately, 
many CAD users return to bench methods. 


The purpose of this workshop is to introduce basic modeling and 
characterization techniques for RF and microwave CAD. Experts will 
discuss how to use simulation to complement the design process for 
wireless communication circuits. The emphasis will be on learning when to 
use simulation to an advantage, how to interpret results, and how to 
establish confidence in the models. Attendees of the workshop will learn to 
trust simulation in order to make engineering decisions. 


A wide variety of topics will be covered, including simulation and modeling 
of passive lumped and distributed components, frequency-conversion 
circuits, and power amplifier circuits. Large-signal device modeling and 
- characterization techniques will also be covered. Identification of the 
appropriate simulation method and model will be covered for each topic, 
and measured versus simulated examples will be given. 


1998 IEEE KADIO AND WIKELESS CONFERENCE (KAWCON’98) 
SUNDAY WORKSHOP 


Modeling and Simulation of Devices and Circuits 
for Wireless Communications 


August 9, 1998, 8:00 am - 5:00 pm 


7 Sana am Registration 
}8:00am | 00 am Introduction 


8:15 am Jeffrey NIST Microwave Network Theory and VNA 
Jargon Calibration 

M/A-COM Stability Analysis of Linear Circuits 
Struble 


10: }10:00 am _| am Break 
10: 15 am | James Sonnet Detailed Analysis of a High 
Rautio Software, Inc. | Temperature Superconducting Filter 
11:15am_ | Stephen Nonlinear Broadband Planar Monolithic 
Maas Technologies, Balanced Mixers and Frequency 
Inc. Multipliers 


eee 15 


1:15 pm Joe Motorola Behavioral Analysis Methods Applied 
Staudinger |Semiconductor | to the Design and Simulation of 
Products Sector | Linear Power Amplifiers 


2:15 pm David Hewlett- Elements of Measurement-Based 
Root Packard Large-Signal Device Modeling 
Company 


55 pm 
3:30 pm Thomas University High-Frequency CAD - The Key to 
Brazil College Dublin | Success in Wireless Design 


4:30 pm Question and Answer Session 


: : 
ke 


pe Che 


Microwave Network Theory 
and VNA Calibration 


Jeffrey A. Jargon 


RF Technology Division 
National Institute of Standards and Technology 
Boulder, CO 80303 USA 
TEL: 303-497-3596 | FAX: 303-497-3970 | E-Mail: jjargon@nist.gov 


RAWCON ‘98 Sunday Workshop 
August 9, 1998 
Modeling and Simulation of Devices and Circuits 


for Wireless Communication Systems 


Outline 


Introductory Remarks 
¢ Microwave Circuit Theory 


Network Parameters 
e Vector Network Analyzers 


Concluding Remarks 


References 


Introductory Remarks 


¢ Definition of a Linear Network. 


¢ Role of Measurements in CAD. 


¢ Importance of Vector Measurements. 


Linear Networks 


Linear Network 
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Role of Measurements in CAD 


¢ Electrical CAD is critical to success. 


¢ Measurements are important in two 
ways: 


Validate elements that can be accurately 
modeled. | 

Provide database for elements that cannot 
be accurately modeled. 


The Importance of 
Vector Measurements 


¢ Complete Characterization of Linear 
Networks 


~* Complex Impedance 
e Accurate Calibration Procedures 
¢ Transformation to Time Domain 
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Microwave Circuit Theory 


Circuit Theory Assumptions 
Low-Frequency Circuit Theory 
Conventional Microwave Circuit Theory 
General Microwave Circuit Theory 
Fundamental Parameters 

Transmission Line Equivalent Circuit 


e The circuit is linear. 


e The circuit is an interconnection of n-port 
components which interact only at the ports. 


¢ The circuit behavior can be predicted solely 
from knowledge of these parameters. 


Two Forms of Circuit Theory 


Low-Frequency 
e Applies to circuits much smaller than a wavelength. 
¢ E-field line integral is virtually independent of path, thus voltage is well-defined. 
¢ Fundamental quantities are voltage and current. 
e Impedance parameters relate v and 1. 


Microwave 


° Circuits may be much larger than a wavelength. 

* E-field line integral depends on path, thus voltage is not well-defined. 
¢ Uniform transmission lines are assumed. 

¢ Fundamental quantities are modal wave amplitudes. 


¢ Scattering parameters are primary, while impedance parameters are defined in 
terms of S-parameters. 
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Conventional Microwave 
Circuit Theory 


Assumptions: 


¢ Mode fields et and ht are real. 
e Wave impedance 71 is constant. 
¢ Characteristic impedance Zo is arbitrary but real. 


Assumptions based on: 


¢ Lossless conductors (loss included as perturbation). 
e TE, TM, or TEM modes. 


General Microwave 
Circuit Theory 


Inhomogeneous Transmission Line Cross Section 


¢ Hybrid modes (not TE, TM, or TEM). 
¢ Common in planar structures. 


Lossy Conductors 


e Skin depth may be large compared to metals. 
¢ Problem at low frequencies. 
e Common in planar structures. 


Complex Characteristic Impedance 


Fundamental Parameters 


Electrical Parameters: 


¢ y: Propagation Constant - signal distortion, loss, and 
delay. 
* Zo: Characteristic Impedance - signal reflection, circuit 
compatibility. 
Physical Parameters: 


* C, G: capacitance and conductance per unit length - 
material parameters of dielectric substrate. 
* L,R : inductance and resistance per unit length - material 


parameters of conductors. 
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Transmission Line 
Equivalent Circuit 


u(Z) LAz RN U(z+A Zz) 


Y= V(GO@L+R\j@C+G)  Zo=./(jal +R) joC+G) 


Characteristic Impedance 
of a CPW Line 


mmmeese REA 2 oy 
Imag{Z ,} 


Frequency (GHz) 


Resistance and Inductance 
of a CPW Line 


Frequency (GHz) 
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Network Parameters 


. © Scattering Parameters 


¢ Traveling Waves 
e Pseudo-Waves 
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Using S-Parameters 


e At high frequencies, S-parameters are measured since it is 
difficult to measure total voltage or current. 


e S-parameters are ratios of traveling waves. 
e They are relatively easy to measure. 


¢ They relate to common measurements such as reflection 
coefficient, insertion loss, and gain. 


e S-parameters can be cascaded for multiple devices to 
predict system performance. 


e They can be converted into Y and Z parameters. 
¢ They are convenient for use with CAD programs. 


Two-Port S-Parameters 
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S-Parameters in Matrix Form 
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Traveling waves ao and bo: 
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Power in terms of traveling waves: 
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Pseudo-Waves 


Pseudo-waves a and b are mathematical artifacts: 
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Scattering and 
Pseudo-Scattering Matrices 
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e Architecture of Vector Network 
Analyzers. 


¢ Calibration Function and Procedure. 
¢ On-Wafer Calibration Issues. 

¢ Verification Techniques. 

¢ Time Domain Network Analysis. 


VNA Block Diagram 


Incident Transmitted 
®) Se 


Source Reflected 
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VNA Calibration: Function 


e Accounts for serious imperfections in the network analyzer: 
Directivity and crosstalk errors due to signal leakage. 
Impedance mismatches in the source, receiver, and load. 
Frequency response errors in the source and receiver. 
Losses in the cables and connectors. 


e Does not account for: 
System drift. 
Repeatability in the switches and connectors. 
Instrument noise. 
Errors in calibration standards. 
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«xem Calibrated Short Circuit 
Uncalibrated Short Circuit 


0 10 20 30 40 


Frequency (GHz) 


Procedure: 


e Measure known standards. 
¢ Process data to determine error coefficients. 


e Correct measured data for device under test. 


Methods: 


¢ Lumped-element standards (OSLT, LRM). 
e Transmission line standards (TRL, Multiline TRL). 
e Electronic transfer standards. 


VNA Calibration: Procedure 
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Forward EXr 


S21 


EDF Slip 22 ELF 
S12 
bo 
EDF = Directivity ETF = Transmission Tracking 
ESF = Port 1 Match ELF = Port 2 Match 


ERF = Reflection Tracking EXF = Leakage 


VNA 12-Term Error Model 


Reverse al ete te ERR 


S21 
S11 §22 


bo a3 
EXR 
EDR = Directivity ETR = Transmission Tracking 
ELR = Port 1 Match ESR = Port 2 Match 


ERR = Reflection Tracking EXR = Leakage “ 


OSLT Calibration 


Most popular form of 
VNA calibration. 


Primarily used in coaxial 
environments. 


Uses lumped-elements. 
Originally intended for 
three-sampler VNAs. 
Variations are possible 
for four-samplers VNAs. 


35 


TRL Calibration 


Next most common form 
of VNA calibration. 


Primarily used in non- Thru | 
coaxial environments. 


Uses transmission lines Reflect _—=_ — 


as standards. 


Requires a four-sampler Line }-—__-—. 


VNA architecture. 


Variations are possible 
for three-sampler VNAs. 


On-Wafer Calibration Issues 


* Known load impedances are rare on planar structures. 

| ¢ Known reflection coefficients are rare, except for 

| transmission lines which have no reflection. 

¢ On-wafer calibration: 
Custom calibration standards correspond to DUT. 
Calibrate out effect of transition and contact. 
Accurately set reference plane and reference impedance. 


¢ Off-wafer (probe-tip) calibration: 


Generic, commercial calibration standards. 
Ignore effect of transition and contact. 


Reference plane confined to the probe tip. 
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Multilime TRL Calibration 


RMR 


Thru 


Multiple Lines 


Reflect 
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Multiline TRL Calibration 


¢ Calibration standards are transmission lines. 

¢ Standards are custom, and also large and numerous. 
¢ Broadband. 

¢ Accurately determines S-parameters. 


¢ Accurately determines propagation constant, which allows 
the reference plane to be set and moved. 


¢ Reference impedance is equal to the characteristic 
impedance of the line. 


Can’t measure pseudo-scattering parameters 
(rer 002): 
Can’t measure impedance parameters. 


Measuring Characteristic 
Impedance 


Often, jaC+G can be estimated. 


When the dielectric loss is very low, thenja@C+G =j@C Bi: 
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Full Electrical Characterization 


RES 


Pseudo-Scattering 
Parameters (50 Q) 


S-Parameters 


Impedance 
Parameters 


Measure Measure 
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OSLT TRL 
postr | ran TRL 


Easy to Implement 


Commercially Available no 
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Resetable Reference yes 
Impedance and Plane 

Complete Transmission 

Line Characterization yes 


Verification Techniques 


Calibration Comparison 


Two similar or different calibrations are performed on 
one VNA, and the error coefficients are compared to 
estimate the magnitude of the maximum differences 
in the S-parameters. 


Network Analyzer Comparison 


Verification devices are measured on two calibrated 
VNAs, and then the S-parameters are compared 
individually for each device, or the S-parameters are 
compared for an ensemble of devices to obtain a worst 
case bound. 
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=_—— Worst-case Bounds 
~——— Measured Deviations 
——— Sliding Load Repeatability 


6 8 10 12 14 16 18 
Frequency (GHz) 


Worst-case error bounds and measured deviations for S11 
between an OSLT calibration and a sliding load calibration. 4, 


VNA Comparison Example 


om Sum of Repeatability Bounds 
~~~ Difference Bound (2 Devices) 
~~ Difference Bound (4 Devices) 


Frequency (GHz) 


Comparison of a three-sampler VNA to a four-sampler VNA. 
Both VNAs calibrated using OSLT with identical standards. 45 


Time Domain Network Analysis 


¢ TDNA is the acquisition of frequency-domain network 
parameters using a digital sampling oscilloscope with time- 
domain reflection/transmission (TDR/T) capabilities. 


¢ This is done by calibrating Fourier-transformed TDR/T 
waveforms using conventional VNA error models and 
calibration methods. 


¢ The significance of TDNA is that it allows laboratories that 
already own lower priced oscilloscopes to make accurate 
frequency-domain measurements. 

¢ The accuracy of TDNA measurements depend on the 
oscilloscope’s settings, drift, jitter, time base distortion, 
voltage distortion, etc. he 


TDNA Measurement Example 


Tota MN M 


10 
Frequency (GHz) 
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Linear measurements are difficult especially 
when the following are present: 


¢ Hybrid modes 

¢ Lossy conductors 

¢ Lossy substrates 

¢ Complex characteristic impedance 
¢ More than two ports 
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The Nonlinear Challenge 


¢ To characterize a nonlinear device, one usually measures 
one or more of the following parameters: load pull, 
spectral regrowth, third order intercept point, adjacent 
channel power ratio, and 1 dB compression point. 


e All of these parameters come from the same nonlinear 
behavior, although they are difficult to compare and 
correlate. 


¢ Our goal: A general nonlinear network analyzer for 
measuring weakly nonlinear devices. 
50 
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Behavioral Analysis Methods Applied To 
The Design & Simulation of Linear Power Amplifiers 


J. Staudinger 
Motorola, Semiconductor Products Sector 
2100 E. Elliot Road, MS EL-712, Tempe, AZ 85284 


Abstract 

This work examines some of the issues and limitations 
of applying behavioral analysis methods to the design 
and simulation of linear power amplifiers driven by 
digitally modulated signals, such as for example, those 
applied in digital cellular handsets. Compared to other 
analysis techniques, this method offers the advantages 
of computational efficiency and a high degree of 
accuracy when certain criteria are satisfied. Such 
criteria are examined and illustrations are included 
comparing the results obtained from this method to 
measured amplifier data as well as to results obtained 
by other simulation techniques. 


Introduction 

Generally, the usual RF design/simulation tools such as 
SPICE and conventional harmonic balance are poorly 
suited for predicting the non-linear performance of 
power amplifiers driven by digitally modulated stimuli. 
In second generation digital cellular systems, the 
modulation bandwidth is small in comparison to the 
carrier frequency. For SPICE based analysis, time 
samples are taken at a fraction of the period of the 
carrier signal which then requires an extraordinarily 
high number of time samples (and simulation time) to 
represent even a short data sequence. The difficulty 
with conventional harmonic balance lies with high 
numbers of frequency components needed to describe 
such a waveform. Several alternatives to SPICE and 
conventional harmonic balance can be applied to this 
problem. One alternative is a behavioral analysis 
approach which is the topic of this paper [1-4]. 
Numerous behavioral models have been suggested and 
a substantial amount of work has been published in this 
area. The focus of this work, however, is to review a 
number of basic issues which affect the success of this 
method when it is applied to predict key power 
amplifier performance parameters, such as adjacent 
channel power ratio (ACPR), output power, and 
efficiency. 


Behavioral Analysis 

To begin, the bandwidth of the modulated carrier is 
assumed narrow (with respect to the carrier frequency) 
and as such the signal can be written as: 


Vin(t) = V (t)cos(@o t+O(t)) (1) 


or alternatively in an I/Q format as 


vin(t) = Re|(Vi(t) + jVo(t) Je" | (2) 
where V(t) is the modulation envelope, 

@, 1s the carrier, and 

©(t) describes phase characteristics. 


Assuming the amplifier non-linearity can be represented 
as a memoryless and bandpass non-linearity, and 
restricting the output voltage to the first spectral zone, 
the output voltage takes the form: 


Vout) =Re| GV(a)e/ CENCE e } 3) 


Where G(V) and (V) describe amplitude and phase 
characteristics of the amplifier as a function of input 
envelope voltage (V). Alternatively, the voltage can be 
expressed in an I/Q format which takes the following 
form: 


NoiyeoRe ee iVa(e)) ‘| joo | 4 
JGAV (t)\(Vi(t) + JVq(t)) 


Where G, and G, are in- and quadrature-phase 
amplitude (only) gain transfer functions expressed as 


Gi(V) = G(V )cos(g(V)) (5) 
and 
GV) = G(V)sin(g(V)) (6) 


In the above equation which describes the output 
voltage around the first spectral zone, the low pass 
nature of the signal is evident. However, it is important 
to emphasize two underlying assumptions, 1) the non- 
linearity & stimuli/circuit must be bandpass and 2) the 
non-linearity must be memoryless. The first requires 
that the frequency-dependence of the amplitude and 
phase characteristics be negligible over several x the 
modulation bandwidth centered around the carrier and 
harmonic frequencies. To a large extent, this is often 
the case for power amplifiers designed for second 
generation cellular systems since the modulation 
bandwidth is small in comparison to the carrier 
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frequency. In addition, circuit aspects of the amplifier, 
such as input & output matching networks applied 
around the device, generally exhibit a wide band 
response (relative to the modulation frequency) thereby 
providing a near constant impedance to the device over 
several x the modulation bandwidth. One exception is 
possible sub-harmonic frequency termination effects 
occurring near the zero’th spectral zone, such as those 
due to drain bias decoupling networks [6,7]. The 
second assumption dealing with memory is often 
violated to some extent due to thermal effects and other 
device memory behavior (e.g., time constants due to 
surface traps). 


Previous work has shown that the behavior model can 
also be utilized to predict both supply current Idd(t) 
(from bias supply Vpp) and power added efficiency [7]. 
Similar to Eq. 4, supply current can be expressed as a 
function of the magnitude ct the input envelope voltage. 
The time averaged value of Iddit) along with the known 
supply voltage allows computing the de power supplied 
to the amplifier (Ppc). 


In general, the amplitude and phase transfer functions 
describing the non-linearity «1G; Vs and 0: V)) can be 
obtained either by measurement ©: simulation (such as 
harmonic balance or SPICE usiny a single tone CW 
stimulus. In the latter case. the active device must be 
represented with an appropriate lurse signa! model, and 
passive circuitry such as the input and output matching 
networks, must be adequately described over a broad 
frequency range to obtain an accurate determination of 
the non-linear gain/phase transte: tunction. 


Issue 1 - Non Linear Transfer Functions 

In general, the transfer functions G and G, are often 
represented as a convergent series ‘odd functions) of 
sufficient order to accurately represent the measured 
am-am and am-pm response ‘The coefficients of the 
series expansion can easily be computed using least 
square techniques. 


An example illustrating the am-am & am-pm response 
of a power amplifier obtained via harmonic methods 
based on a single tone CW stimulus 1s shown in Figure 
1. Consider representing G, as a series expansion of 
Bessel and Polynomial functions taking the form: 


2M +1 

G,= Ya, JUkAE) (7) 
k=odd 

and 
2M +1 

G,= > b, (AE) (8) 


k=odd 


where A is the amplitude of the signal envelope and & is 
a normalizing function. The residuals from fitting G; to 
the am-am & am-pm data for values of M ranging from 
2 to 40 is shown in Figure 2. A Q/R matrix factorization 
method based on a least squares method was applied to 
extract series coefficients. Clearly, as the order of the 
series increases, the residuals decrease and the series 
expansion more closely represents the data. An 
exception is noted for the Polynomial series at high 
orders (M>20) where the residuals begin to increase 
due to high condition of the matrix. It is also noted that 
for a given value of M, the residuals are generally 
smaller for the Bessel series. Although not shown, the 
results for the quadrature gain transfer function (Gg) 
exhibit similar behavior. 
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Figure 2 Residuals from fitting G; to series expansions 
of Bessel and Polynomial functions. 
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A comparison of the amplifier’s IM3 performance 
computed using harmonic balance contrasted to that 
obtained by behavioral analysis methods is illustrated in 
Figure 3. An inspection of Figure 3 suggest that in 
applying behavioral analysis techniques, the results 
obtained with transfer functions G; and G, expressed 
with Bessel functions result in more accurate 
predictions of IM3 compared to Polynomial expansions. 
It is also noted from an inspection of Figure 2, at M 
values of 10 and 15 the former series expansion has 
significantly smaller residuals than the latter. The 
results also suggest that for increasing values of M, the 
behavioral model predictions become more accurate. 
This trend is in agreement with previously reported 
results [6,7]. 


Issue 2 - Sub-harmonic Terminations 

Previous results [5,6] have shown that sub-harmonic 
frequency terminations on the active device’s gate & 
drain terminals on or around the frequency of the 
modulation envelope can significantly effect linearity 
performance of the amplifier. This effect is illustrated 
by way of example. Consider a 1 W power amplifier as 
illustrated by the partial schematic shown in Figure 4 
where a 250 nH inductor is utilized for drain biasing. 
The IM response of the amplifier when driven by a two 
tone sinusoidal signal with tone spacing of 10 kHz and 
1 MHz is calculated using harmonic balance methods 
and is illustrated in Figure 5. In addition, the amplifier 
am-am & am-pm response is calculated using harmonic 
balance by sweeping the power of a single tone CW 
stimulus. Both in-phase and quadrature-phase transfer 
functions are then calculated and the behavioral analysis 
method is exercised to determine the amplifier’s IM3 
performance for the above mentioned two tone 
sinusoidal stimulus at the two tone spacings. These 
results are also shown in Figure 5. The amplifier’s IM3 
performance computed using the two analysis methods 
are nearly identical for the 10 KHz tone spacing. 
However, at the higher tone spacing of 1 MHz, the two 
analysis methods differ substantially over the entire 
power range. It is also observed that the upper and 
lower frequency IM frequency products (including 3", 
5", and 7” ) are unequal in amplitude. It is also 
observed, that the behavioral analysis results do not 
change with tone spacing. 


For this case, it will be shown that these differences are 
due to device sub-harmonic termination effects related 
to the 250 nH drain bias inductor. At envelope 
frequencies of 10 kHz and 1 MHz, the 250 nH inductor 
presents an impedance at the drain terminal of .016 and 
1.6 ohms, respectively — both of which are a relatively 
smal! value. 
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Figure 3. Amplifier IM3 response computed using 
harmonic balance (HB) compared to 
behavioral model predictions a) w/ G, and 
G, expressed as series expansion of Bessel 
functions, and b) w/ G; and G, expressed as 
series expansion of Polynomials. 


The amplifier is now modified as shown in Figure 6 by 
including series resonant circuits connected at the 
device’s intrinsic drain terminal. The resonant 
frequency of these circuits are set to coincide with the 
fundamental and harmonics of the envelope frequency, 
AUN Ieavelopen N= be, oe -2 8 (€:g., 1; 2,°3 2. MHz). It is 
noted that these resonant circuits have no influence at 
the carrier frequency and harmonics of it. The 
amplifier’s single tone am-am and am-pm response 
remains unchanged. The IM _ performance is 
recalculated using harmonic balance methods and the 
results are found to be nearly identical to those shown in 
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Figure 5, matching the results for the 10 KHz tone 
spacing and agreeing very well with behavioral analysis 
predictions. 


250 nH 


Figure 4 Partial schematic of a 1 W power amplifier 
where a 250 nH inductor is utilized for drain 
bias de-coupling. 


Figure 6 The amplifier is modified by including series 
resonant circuits connected at the device’s 
drain terminal with the resonant frequency 
set to coincide with the fundamental and 
harmonics of the envelope frequency thereby 
presenting a zero ohm termination to these 
frequencies components. 


This effect is investigated further by examining the 
envelope signal at the outputs of both amplifiers around 
the first spectral zone (Figure 7). The envelope signal 
corresponding to amplifier using a 250 nH drain bias 
inductor (Fig. 4) exhibits an envelope signal which is 
asymmetrical wrt to time with the asymmetry due to 
envelope signal contributions. The envelope signal has 
the appearance of a memory effect even though the 
large signal model utilized in the simulation is ideal in 
that it is memoryless (i.e., the instantaneous drain 


current does not depend on prior gate-source and drain- 
source voltages) The other amplifier configuration 
(Fig. 6) exhibits an envelope signal which is highly 
symmetrical wrt to time. 
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Figure 5 The IM3 performance of the power amplifier 
is simulated using harmonic balance (HB) 
and behavioral analysis methods for an input 
two-tone sinusoidal signal with tone 
separations of 10 KHz and 1 MHz. a) 10 
KHz tone spacing, and b) 1 MHz tone 
separation. 
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Figure 7 The envelope signal at the amplifier’s output 
terminal around the first spectral zone. a) 
Envelope signal corresponding to amplifier 
configuration of Fig. 4, and b) Envelope 
signal corresponding to amplifier 
configuration of Fig. 6. 


Summary 

Behavioral analysis methods can accurately predict 
distortion characteristics, output power, and efficiency 
of power amplifiers driven by multi-tone sinusoids and 
digitally modulated signals when certain underlying 
criteria are satisfied. In general, model accuracy 
improves by representing the gain transfer functions 
with series expansions of increasing order. It is shown 
that sub-harmonic frequency terminations, such as those 
related to drain bias networks, can significantly effect 
the linearity performance of the amplifier. Such effects 
are not modeled by the behavioral analysis methods 
presented in this work. 
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Summary 


e Overview and Context 
¢ Non-linear Device Modelling 
— ¢ Circuit Simulation Techniques 
e Systems-Level Simulation 
e Future Challenges 
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UMTS 
(Universal Mobile Telecommunications System) 
IMT- 2000 
(International Mobile Telecommunications - 2000) 


e UMTS is the term commonly used within Europe to describe the 
next generation of wireless, wideband digital communication 
systems with enhanced characteristics such as full-motion video, 
graphics-intensive Internet, flexible service provision, universal 
coverage, and support for data rates up to 144kb/sec (wide-area) or 
up to 2Mb/sec (indoor/local); 

¢ Ata global level, UMTS represents Europe’s input to the effort by 
ITU-R to standardise third generation wireless systems, known as 
IMT-2000 , previously referred to as FPLMTS. The ‘2000’ signifies 
both the approximate frequency band (2000MHz or 2GHz), as well as 
the approximate initial date for deployment. 


University College Dublin 


Evolution Towards UMTS 
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Standardisation of Air Interface for UMTS 


In January 1998, ETSI agreed a critical milestone for the future 
development of UMTS, namely, definition of a single standard for 
the air interface, hoping to allow Europe to build on its GSM 
success; 


The chosen concept uses WCDMA for wide-area applications and 
TD/CDMA for low-mobility high-capacity indoor applications, and 
will be jointly promoted by leading telecommunications suppliers, 
including significant support from Japan; 


The WCDMA decision was strongly supported by operators who will 
not need new cites or cell structures - they can simply overlay 
existing GSM networks in migrating to the new WCDMA hardware, 
which will operate to open standards; 


The standard will support basic speech as well as flexible high 
speed data (circuit-switched - e.g. ISDN and packet e.g. TCP/IP), with 
Ressible multiple simultaneous services. 
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some Technology Features of UMTS 


¢ Use of adaptive antennas (on MS and BS); 

¢ Integrated satellite component (S-UMTS) for very-wide 
area operation on a global scale, ultimately maybe linked 
to GPS; 

e Extensive use of IN concepts 

¢ ‘Software radios’ may permit highly flexible terminals, self 
customising to local standards or capable of down-loading 
upgrades, new features etc. 
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Emerging MM-Wave Markets 


e The characteristics of mm-wave which make it problematic for 
traditional microwave relay applications, actually make it ideal 
for short distance links, using cellular-type concepts; 

e Multiple frequency re-use becomes possible over very short 
distances (e.g. indoor broadband wireless, pico-cell concept at 
60GHz, taking advantage of strong absorption in this range); 

e New volume markets are under development in LMDS around 
29GHz, MVDS (42GHz2), broadband wireless LAN’s, satellite links, 
MBS (60GHz), automobile ACC-CW/CA (77GHz) etc. 
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From Cordless to Mobile Broadband (MBS) 


Wired Terminal 
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Device/MMIC Technologies 
Silicon: 
¢ bipolar junction transistor (BUT); 
¢ MOS-based (LDMOS, RF CMOS..); 
Silicon Alloys: 
¢ SiGe Heterojunction Bipolar Transistor (HBT). 


Compound Semiconductors: 

e GaAs FET (MESFET); 

¢ High Electron Mobility Transistor (HFET/HEMT); 
e Pseudomorphic HEMT (PHEMT); 

e GaAs-based HBT’s; 

e InP-based FET’s. 
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Block Diagram, Typical Multi- 
Technology RF Transceiver 
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Technology Speed Figure of Merit 


vs Time 
ft Hemts,HBTs —_, -_. 
a5 : oe . 
100GHz| ie Pd 
30GHz et et, 0 25u 
~ oy . too ee ne ere tenn rier eee core 


3GHz| GaAs pric ne a 


we ee re we) 7. 0.5u 
| 10GHz Ses . 7. a 0.6u 

| eae (oe ‘ 0.8u 

| 

| 

| 

| 


Bipolar a, 
MOS 
ef 
| | | | | | 
75 nA 73 81 83 85 87 89 31 


93 95 97 39 
= oS 
Year | 


University College Dublin 


GaAs-Based Technology Trends 


Breakdown of MMIC Designs by Technology 


1993-1995 2000 


Figure 2 
The MDE Consortium's view of the increasing usage of PHEMT and HBT circuits in 


MMICs, indicating that they will likely dominate the field by the year 2000. In 
response, the MDE project is focusing on improving microwave and millimeter wave 
non-linear region modeling. 
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High-Frequency CAD 


e CAD is now a critical enabling technology across almost every aspect of 
microwave and RF applications, but especially in MMIC design; 


e CAD may be applied at many different levels, e.g.: 


Process Simulation: predicting distributions of dopants, layer 
structures... from given growth conditions; 

Physics-Based Device Modelling: predicting/optimising performance 
(and ideally yield) from doping, geometry and materials information; 
Electromagnetic Solvers: analysing or even synthesising complex 
physical geometries with Maxwell's equations; 

Circuit-Level CAD: Analysis/design of linear/non-linear microwave 
sub-systems for gain, noise, efficiency, linearity, output power etc.; 
Systems-Level CAD: predicting the performance of complete 
microwave systems, possibly studying trade-offs between 
baseband/DSP solutions and RF techniques; 
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General Trends in High Frequency CAD 


shift from military to commercial markets; 
shift from Workstation to PC; 


increased complexity of CAD tools/interfaces - role of vendors in 
contract design; 


high cost of ownership of mainstream CAD environments but 
these are then often ‘under-used’; 


CAD a strategic technology (e.g. for UMTS/IMT-2000); 
Internet as a means for sales, marketing and support; 
links between high-frequency, DSP, optoelectronics.. 
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some Technical Trends in High-Frequency 
CAD 


¢ increasingly stringent system specs. placing severe demands on 
existing CAD tools; 

¢ volume exploitation of mm-wave markets (LMDS, MVDS ..) 
requiring better mm-wave models, yield analysis tools, 
EM/physics-based simulation...; 

e users need more synthesis tools as well as tools for analysis: 
increasing interest in ‘design assistant’ software; 

e sub-micron CMOS an increasingly competitive RF technology at 
the lower frequencies; 
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The European EDGE Project 


e EDGE is a project part-funded under the European Commission 
ESPRIT programme, in the area of high-frequency analog CAD; 


¢ It is about 2 years in duration and due to finish during 1998; 


¢ The total buget for the project is 2.8 Mecu of which the EC 
contribution is 1.6 Mecu, representing a total work effort of about 
27.5 Person-Years; 


¢ The project consortium consists of a mixture of software 
developers/vendors (5) and users/component manufacturers (3). 
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The EDGE Consortium 


¢ Co-ordinator: University College Dublin (UCD - IRL) 
|. (National University of Ireland, Dublin) 
e Partners: Barnard Microsystems Ltd. (BML - UK) 
Dassault Electronique (DE - F) 
GaAsCode Ltd (GC - UK) 
GEC-Marconi Materials Technology (GMMT - UK) 
Jansen Microwave GmbH (JMG - D) 
Philips Microwave Limeil (PML - F) 
¢ Ass. Partner: Universita di Roma Il ‘Tor Vergata’ (UR - /) 


[vendor/developer] [user/manufacturer] 
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Specific Objectives of EDGE Project 


e Provide Improved Working Links between Commercial European CAD 
Tools; 


¢ Provide Useful Links between European Tools/Models and US-vendor 
software; 

¢ Support Advanced User Interface Development for MMIC Design and 
Evolution of a Common User Interface Style; | 

¢ Integration of Time-Domain, EM analysis and Yield Enhancement CAD 
tools for MMIC Design; 

¢ Provision of a Standardised Advanced Non-linear European FET 
Foundry Model (the COBRA model); 

¢ Develop a future collaborative business structure for commercial 


exploitation of European CAD expertise both within EDGE and outside 
(e.g. from universities/institutes). 
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Three Main Axes of Technical Development 


e Semiconductor device characterisation, modelling 
and parameter extraction; 


e Passive component modelling and electromagnetic 
simulation; 


e Simulation/optimisation engines and solution 
strategies and algorithms. 


‘from Gate recess depth to Bit Error Rate’ 
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¢ Nonlinear Device Modelling 


* Physics-based non-linear models; 
* Equivalent Circuit Based device models; 
* Table-Based Models, other approaches ... 
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Physics-Oriented Nonlinear Modelling 


In this approach, the physical behaviour of the non-linear device is 
encapsulated in a set of equations, which assume a knowledge of 
basic material parameters for the device material, and full 
geometrical! and process information for the device under study; 
The potential advantages are clear:- clear physical insight into 
behaviour, potential for process optimisation for performance and 
improved yield; no need for fitting large numbers of parameters to 
empirical models; automatic consistency over operating regions 
and predictive capability; 

Advances in computing power at low cost are bringing this kind of 
simulation more-and-more within the reach of circuit designers; 
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Hydrodynamic Equations 


the standard drift-diffusion equations of semiconductor physics are 
inadequate to describe hot-electron transport and velocity 
overshoot in short-dimensional structures; 

carrier transit times in such structures are of the same order as the 
energy and momentum relaxation times, and the standard local 
models of mobility break down; 

in such cases, a more rigorous solution is required based on the 
Boltzmann transport equation, applying conservation of particle 
number, momentum and energy, as the zeroth-, first- and second- 
order moments of this equation, respectively; 
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Basic formulation of hydrodynamic equations: 
(n = electron concentration; p= momentum; W= energy density; 
u = heat flux; subscript ‘c’ indicates collision contributions) 


any -(nv) = (=) 
at at). 


Pe VW -p+p-VV= -gnE-v(nT)+(2) 


(VW) = -gnv-E-v.(vat)-v-0+(SF) 
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Constitutive Relations (parabolic energy bands): 
(k = thermal conductivity of electron gas) 


2 
U =e = pein eV nT + 2 njv] 


Heat Flow: 


(T,, = lattice temperature; ¢ = specific heat; p = density; x = lattice thermal conductivity) 


OT, 
ot 


cp =V.(x,VT)+.6+4E,( 2] 


ot 
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some Issues in Hydrodynamic Analysis 


e There are many possible variations in the physical equations, 
depending on effects considered or assumptions made (cf. 
electromagnetic analysis where the equations to be solved, Maxwell’s 
Equations, are much better defined); 

¢ Non-parabolicity of the energy bands is probably necessary but 
introduces much additional complexity and uncertainty; 

e In devices involving quantum wells etc., the equations can be extended 
to include modes arising from solutions to Schrodinger equation at the 
cost of extra complexity (‘quantum hydrodynamics’) 

e Determining values the basic physical parameters may be problematic 
over the wide range of parameters required and ordinary users do not 
usually have access to detailed process data; 

¢ Numerical analysis of such equations under large drive conditions is a 
non-trivial problem and still time-consuming 
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Non-Linear Equivalent Circuit Modeliing 


This is generally the most popular and well-understood approach among 
designers; 


it is customary to divide the model into an extrinsic part, representing the 
electrical (and possibly thermal) environment between the active part of the 
semiconductor and the accessible outside terminals, and an intrinsic part, 
describing the core device behaviour; 


While broadly justifiable on the basis of physics, the usual intrinsic model 
formats represent a major simplification of what is a much more complex 
underlying behaviour at microwave frequencies - e.g. the conducting 
channel may behave in reality more as a distributed, active R-C line, device 
fingers are themselves distributed lines at high frequencies etc. ; 


The extrinsic network is dominated by parasitic elements, and may become 
complicated with element values which are hard to determine, especially 
over broad operating bandwidths (~ f,) and if the device is packaged; 
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Topology and Parasitics 


Intrinsic 
Transistor {| 
Model | 


accurate extraction of parasitics is critical! 
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Some Issues in Non-linear FET Modelling 


Continuity and high-derivative accuracy of model equations; 


‘Charge-conservation and non-quasi-static behaviour of 
Capacitance; 


Dispersion / non-quasi-static behaviour of drain current / 
bias-dependent delay/ current lag effects; 


Scalability up to very large device sizes; 


Subthreshold, avalanche breakdown, behaviour around 
Origin and into third quadrant; 


Thermal effects and self-heating; 

Distributed behaviour at very high frequencies; 
Parameter extraction and measurements; 
Implementation in CAD environments; 

Model Validation. 
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Calculated derivatives of the drain current model function 
(0.2x120Um PHEMT: Vpo = -1.0V; Vds = 3.0V) 


COBRA Model _ Modif, Materka Model _ 
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Calculated derivatives of the gate-capacitance model function 
(0.2x120Um PHEMT: Vpo = -1.0V; Vds = 3.0V) 


_Materka Model 
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DC Characteristics - 0.5x(4x75)um MESFET (GMMT - F20) DC Characteristics - 0.5x(4x150)um MESFET (GMMT - F20) 
(Vgs =- 1.6Vto +0.4V; Vpo =- 1.1V) (Vgs =- 1.6Vto +0.4V; Vpo =- 1.1V/) 


i 


"| © measured] 
. [== modelled | 


woe 4 Vds [V} 


DC Characteristics - 0.5x(6x200)um MESFET (GMMT - F20) 
(Vgs =- 1.6V to +0.4V; Vpo =- 1.1V) 


DC Scaling 
Tests for 
EDGE Non- 
Linear FET/ 
PHEMT 
model 
(COBRA) 


Vds[V] 
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Simple FET Capacitance Model 
Vv 


eum, Gate 


Total static Gate charge (partition line assumed bias-independent): 


Qe (Ves Ven) = Qgs(Ves) + Qgp(Vep) 


This leads to a simple (quasi-static) capacitance model, as used 
in the early Curtice model etc. E.g.: 


Creo Codo 
C,.(Vs) = 3 Cya(V ga) = 
Clic Vgs/o,]" ide Vga/i]™ 
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Charge Conservation 


¢ Experimental results and physical intuition suggest that the simple model is 
inadequate, as C,. and C,4 will depend significantly on both bias values, 
especially below the ‘knee’ or in the ohmic region of the drain 
characteristics; 
¢ Many simulators then require use of quasi-static partitioned charge 
formulations of the form: 
Qgg(Vgs(t), Vas(t)) and Qgp(V,q(t), Vag(t)) 
¢ Each of these is assumed to be a state variable which conserves charge, in 
the sense of returning to the same initial condition after having been cycled 
through an arbitrary closed dynamic path in the V,,., V5 Plane: 
V C2 
gs 


C1 
Q, Vg 
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¢ The equations of Statz et al. [ED Feb. ‘87 p.160] represented an 
important step in improving the description of FET non-linear charge by 
introducing a semi-empirical formulation for the total Gate charge. This 
formulation is not conservative in terms of Q,, and Q,,. However, it may 
be argued on physical grounds that these charges are indeed not 
conserved in FET’s under dynamic operation. The following example is 
quoted by Statz et al., based on a starting charge distribution 
corresponding to Vg, and V,: 


Vg (< OV) 


Gate 


V,(> OV) 3 OV) 


oh ey (2) Vp, next decreased to V,, 


(1) Vg, first increased to V., Net charge transfer 
from *@.= tes On, 
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Simulator Issues for Charge Modelling 


¢ Statz and similar equations create problems in the context of the conventional 
equivalent circuit topology when implemented in charge-based simulators (such as 
SPICE) as the lack of charge conservation causes spurious current components to be 
generated; 


¢ This can be corrected by adding additional controlled current sources or 


‘transcapacitances’ into the model. E.g. given a charge function Q(V,, V,), we can write 
the full differential as: 


dQ=C.dv, +C€.dV, = coed ay, 
OV, OV, 
It is readily proven by vector algebra that ‘conservation’ or path-independence is 
assured if: edie) é 
OV> OV, 


* Note that some simulation strategies (e.g. state-space solvers) can operate directly on 
the nonlinear capacitances C,, and C,,. Strictly speaking, these kinds of solution also 
cannot be conservative, because of the time-discretised nature of the numerical 
algorithm used, although with small time-steps and good algorithms, spurious effects 
of this kind can be kept to a minimum. 
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Dispersion in Drain Current 


The conventional approach to FET nonlinear modelling regards charge 
as quasistatic but has to contend with the well-known discrepancy 
between dynamic and static drain current characteristics (even at 
constant temperature); 


These ‘dispersion’ phenomena are thought to orginate from channel/ 
interface traps and surface charge, and are characterised by rather low- 
frequency relaxation times ( 1kHz - 1MHz). They are process-dependent 
and often more severe for MESFET’s than for HEMT’s; 


The effect is that primarily drain-source conductance is affected, by 
being increased significantly, but so also is transconductance, by being 
reduced; 

Pulsed DC measurements are a valuable tool for exposing this 
behaviour as described in a number of recent papers. However, it can 
be difficult experimentally to separate thermal and ‘dispersion’ 
phenomena; 
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A first-order dispersion correction is 
possible by adding an RC branch in 
parallel across the (quasi-static) 
drain current; 


To correct for g,, and g,, over a wide 
bias range, a second nonlinear 
generator may be added used as 
shown; 


Even better adjustment is possible 
by also adding an extra noniinear 
conductance g,; 

The determination of C, and its bias- 
dependence is facilitated by pulsed 
DC measurements. 
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Thermal Issues in Device Modelling 


pulsed, multi-bias DC/S-parameter measurements are sometimes 
made, since results are then expected to be more appropriate to the 
temperature at the given Q-point 


in reality, bias point (and temperature) change with increasing RF 
power levels, and also as circuit conditions change 


full electro-thermal non-linear simulation is needed properly to 
account for all the interactions involved 
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COBRA Model 


COMET FET Model Extractor — 


reeseemene 


Single-tone large- 
signal test results 
(9 GHz) fora 
0.2x(6x50)um 
PHEMT, using the 
COBRA model 
scaled from 
0.2x(4x30)um 
{Vpo = -1.0V). 


Pout [dBm] 
Pout [dBm] 


‘20.0 15.0 -100 -5.0 0.0 5.0 10.0 
Pin [dBm] 


b) 


Vgs Vds 30.0 
a) -0.9V, +0.6V; 10.0 


b) -0.9V, +6.0V; 2°” oe 


-20.0 


c) -0.6V, +3.0V; -40.0 fom 


Pout [dim] 
Poult [dBm] 


d) 0.0V, +0.6V. =: 


-20.0 -15.0 =10.0 -5.0 0.0 5.0 10.0 
Pin [dBm] 


c d 
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Modelling of HBT Devices 


e The HBT is a relatively complex device, for which there is now 
considerable demand for effective CAD-oriented modelling 
support; 


e GaAs-based HBT’s have been under development for a longer 
time and have illustrated some of the difficulties in adapting 
standard BJT models for linear or nonlinear modelling; 


e SiGe HBT’s show significant differences: the thermal issue is 
important but perhaps not so critical, however base push-out 
(the Kirk Effect) is considerably more important; 

e Users require compact, robust model implementations in 
commercial simulation environments, backed up by clear 
parameter extraction procedures. There is still much work to be 
done to meet these needs. 
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Problems Applying the Gummel-Poon Model to 
lll-V Based HBT’s 


e Many devices do not exhibit any useful ‘constant current gain’ 
region, yet the existence of such a region is a fundamental 
principle of the GP model; 

e In the high-current region of the Gummel Plot, the measured 
curves tend to deviate from linearity, but for a fixed V,, these 
deviations are generally unequal for the HBT, rather than ~equal 
as predicted by GP; 

e For HBT’s, the deviation from linearity for |, may be larger than 
that of |, or the other way round; 

¢ Temperature must be solved for explicitly rather than being 
assigned in advance. ‘Gain collapse’ due to thermal instability 
can occur in even in non-power processes. 


University College Dublin 


AlGaAs/GaAs HBT Characteristics 


Solid lines: simulated results: 


P : 
j ne ‘ aaa Circles: measured results: 
3 a i 
3 Selo” 3 


Leakage Current —— fs af . Res 
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Region 


BE a eg 
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Transistor Nonidealities in Static GP Model 


Recombination 
Leakage 


Base-Width 
Modulation 


(Early Voltage) 


High-Level 
Injection 
(Webster Effect) 


Base Push-Out 
(Kirk Effect) 


Local Thermal 


Effects 
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Gain Advantage SiGe vs. BJT 


AE 


- g.Ge(X=0) 
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Emitter 
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Collector 
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g,Ge 
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Kirk Effect (Base Push-Out) 


. E(x) 
When a bipolar is operated such 


that high-level injection effects Fie Base 
occur in the collector, minority eis 
carriers are injected into the 

collector region adjacent to the 

base, widening the electrical 

width of the base and leading to a 
degradation in current gain and 

high frequency performance; E(x) 


This effect is equivalent to ’quasi- 
saturation’ where the internal 
base-collector metallurgical 
junction becomes forward biased 
while the external base collector- 
terminals are reverse biased; 


Emitter 


Collector 


Collector 


Current-Induced Base 
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Basis of MMIC Sensitivity Analysis and Yield 
Forecasting in EDGE 


Material and 


Device structural 4 


Circuit 
Simulator 


OUTPUTS: 
mic oD 
Sensitivity and | 
Yield 


Frequency 
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Test of Yield Forecasting Software In EDGE 


LOW_NOiSE.AMPLIFIER) Gain (dB) 


Sharp Tail on 
fon Profile 


} Measurement results 
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He First Version of FetWin 
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FetWin Yield Analysis Channeling — 
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Comparison between Physics-Based Yield Predictor 
(FetWin) and measurements of 100 samples: histogram 


of S21 (dB) at 10 GHz. 
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e Simulation Techniques 


* Time-Domain Methods; 

* Volterra Series; 

* Harmonic Balance Techniques; 
* Circuit Envelope; 

* Convolution Methods; 

* Behavioural Models 


Non-Linear High-Frequency Simulation 


The central problem is that the /inear parts of the system are best 
described in the frequency-domain (S-parameters...) whereas the 
non-linear parts tend to produce integro-differential equations in 
the time-domain. The numerical challenge is to find methods 
which optimally handle both domains; 


The many formats adopted for digital wireless communications 
present particular challenges in calculating efficiency, very low 
distortion levels, ‘adjacent channel interference’, etc. mainly 
because of the spectral resolution and range required (e.g. ~ 
90Hz resolution up to many harmonics of 2GHz, perhaps with 
multiple simultaneous tones) 
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TIME-BOMAIN & REAL FREQUENCY-DOMAIN 
RELATIONSHIPS . 
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Time-Domain Methods (SPICE-type) 


Circuit is described using a nodal-numbering scheme, and a wide range 
of ideal linear and non-linear elements are possible, including ideal 
transmission lines; 


DC solution involves formulation of circuit equations using Modified 
Nodal Analysis as a non-linear algebraic system which is solved 
iteratively (LU decomposition, pivoting and sparse matrix techniques). 
Non-linear resistive elements are transformed to a linearised 
(source/resistance) "companion" form at each iteration. 


For transient analysis, the solution interval in the time domain [0,T] is 

divided up into distinct time-steps, t,, t,,,... The Sparse Tableau 

differential equation formulation for the vector x of solution variables, is: 
BUX 0G tO eeX OX 


the dimension of the vector is (2.b+n), where b= number of branches, n= 
number of nodes; 
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This system is solved numerically using a single-step implicit scheme 
(e.g. trapezoidal), starting from the DC solution. This implicit 
formulation effectively reduces the non-linear dynamic circuit functions 
(including charge and flux) to equivalent non-linear resistances. Thus 
transient analysis becomes a sequence of non-linear algebraic 
solutions, or, in effect, repeated DC solutions; 


Problems with "stiff" systems (widely-varying time constants) make it 
desirable to incorporate some form of adaptive time-step control, based 
on meeting a local truncation error tolerance. Unfortunately, the 
stability properties of the algorithm become very difficult to predict in 
this situation; 


If a steady-state solution is of interest, the solution must be run until the 
initial transient period has died away, and then Fourier Analysis may be 
performed over a cycle via the FFT to compute harmonic levels. 
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Advantages of Time-Domain 


Solution starts from DC and follows a physically-consistent trajectory in 
phase-space. If parametric effects, spurious oscillations, or even chaos 
are possible they can appear in the solution; 

Very general excitations in the time-domain are possible; 


Solution is usually quite accurate and the numerical procedures show high 
level of robustness; 


Disadvantages of Time-Domain 


Method can be very inefficient if only a steady-state answer is of interest, 
and a long transient must be solved. Variations such as the "shooting" 
method may help, but they are complex and not widely available; 

The adaptive time-stepping severely compromises the ‘noise-floor’ of the 
FFT - important in sensitive calculations of distortion products; 


Description of realistic, arbitrary distributed external circuit microwave 
environments is not possible. 
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Harmonic Balance 


Introduction 


This method is based on the assumption that the (forced) non-linear 
system has a unique periodic solution, and that both the nature of the 
excitation, and the degree of non-linearity, are such that a relatively small 
number of harmonics (K) can describe all waveforms completely; 


In practical application of HB, an iteration takes place in which a trial 
solution vector (with each element described by a truncated Fourier 
Series), is passed on the one hand through the linear external network in 
the frequency domain, and on the other hand through a non-linear sub- 
system described in the time-domain. The objective of the iteration is to 
adjust the unknown coefficients in order to "balance" all harmonic 
components, so that the final representation is consistent in each domain; 


If a good initial guess is available, very rapid determination of the steady- 
state solution is possible. The method is is offered within all commercial 
software environments, and is very widely used in non-linear CAD. 
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Basic Method of Harmonic Balance: 


Let P(T,) be the set of bounded periodic waveforms in T,. Let ue P(T,) 
be a stimulus waveform in a DE of the form (f continuous and real): 
f(x, X,u)=0 


Assume that the solution xX exists, is real and x ¢ P(T,). Then: 


x(t)= + X(k -e/Kot and u(t)= Saree etvnere” gem ee 
Kk k 10 ali 
=—co =—oo (e) 


Substitute the assumed solution and its derivatives into f, which must 
now also be € P(T,). It can therefore be written as a Fourier Series: 


f(x(t), x(t),u(t))= ¥ F(X,U, k)-e/K°o" 
where: X=, X(- 1), X(0), X(7),-.-) 
U=[--,U(-1),U(0),U(1),--] 


Solve the following system of non-linear algebraic equations for X: 
FIX UKHO V ke Zs 
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Formulating the Harmonic Balance Equation 


Linear 
Sub-Circuit 


Let Y be the linear circuit admittance matrix, /, be a vector of Norton 
excitations, and Q(V) and /,(V) represent non-linearities, then we 
require a solution for V= vector of (unknown) voltages at Ports 1 -- N: 


E(V) = I, + Yuu V + {2.Q(V) + 1,(V) = 0 
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Solving the Harmonic Balance Equation 


A great deal of research effort has been directed towards finding 
efficient, reliable methods for solution of the HB Equation. Possible 
methods include optimisation, splitting, reflection... A conventional 
technique with good convergence properties from a close initial guess, 
is that of Newton or quasi-Newton: 


yer) yO) _ a. Ee . E(v) 


a, is a scalar relaxation parameter ¢ [0, 1], and J() is the Jacobian of E 
with respect to V. However. with several ports and numerous 
harmonics, the matrices involved can become extremely large. Note, 
also, that high accuracy in computing the Jacobian is very desirable 
(analytical formulations are preferable, wherever possible); 
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Krylov Techniques 


e Up to recently, HB methods struggled with the numerical complexity 
required to solve large circuits with many tones as required by digital 
wireless applications; 


e An important recent advance has been the application of special 
methods for solving very large, relatively sparse, asymmetric equations 
of the form A.x = b. In particular, methods have been proposed based 
on generating a basis for the so-called Krylov sub-space, spanned by b, 
A.b, A?.b, --- A*«-).b, One example is the Generalised Minimal Residual 
Method (GMRES - see below); 


¢ These methods have been recenily introduced into commercial 
simulators (ADS - HP-EEsof; Serenade 7.5 - Ansoft). There are reports 
of successful analysis of circuits with 1000’s of tones and over 100 
active devices, each described by a full non-linear equivalent circuit 
model. 
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Generalised Minimal Residual Method (GMRES) 


e Start with an initial iterate x,, using a Preconditioning Matrix H which is 
an approximation to A’: x, =H.b. The residual error is then [b - A.x,]. 


¢ It can be shown that successive iterates can then be constructed using 
H, so that after k iterations, the solution estimate is: 


Xx 2 xX, + H-[r, +(1— AH) -r, +--+ (I- An)" r,| 
with residual: r,, = (I- AH)*-r, 


In GMRES, coefficients a, a,, a%,., are added as follows which are 
adjusted to minimise the Euclidean norm of r,, : 


xX, = Xo+H-|ag-t, + oy (I= AH): ty +o + (4): (I= AH) «4, | 


To achieve this, QR decomposition can be used to construct an upper 
Hessenberg matrix of size (k+1) x k. Using k ~ 10, means the problem is 
of manageable proportions and easy to solve. 
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Sources of Error in Harmonic Balance 


(1) Using Insufficient Harmonics. /t can be difficult to decide a priori what is an 
approriate number for K. Using too many harmonics is inefficient. 
Engineering judgement, or even approximate analytical solutions can give 
guidance and also provide good solution-vector starting estimates; 


(2) Aliasing Errors in the FFT. The FFT is an exact method for computing 
Fourier Series, but only for band-limited signal spectra and subject to the 
Nyquist condition. Practical non-linearities do not preserve band-limiting, 
and then the FFT causes neglected high-frequency components to appear 
as error components on calculated (especially high-valued) harmonics. 
This problem can be reduced by retaining sufficient harmonics and over- 
sampling in the time-domain to preserve accuracy in calculated levels; 


(3) Iteration Closing Error: This arises from the non-exact, numerical solution 
of the HB equation. 
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Advanced Harmonic Balance 
¢ A number of generalisations of Harmonic Balance have been developed 
in recent years, which allow the method to be used in a wide range of 
applications, e.g. mixers and oscillators. Another important case is 
multi-tone excitation of amplifiers, where several strong tones are 
applied simultaneously, at perhaps incommensurate frequencies; 


¢ Where the frequencies are commensurate, the most straightforward 
approach is to use a fundamental period sufficiently long to produce 
periodicity. Very long-record FFT's may be required however; 


¢ If the frequencies are incommensurate, it is generally still possible to 
find a time period over which the waveform is approximately periodic. A 
number of "almost-periodic" transform (APT) techniques have been 
developed for this case. For example, representing the voltage by: 

K-1 

v(t)= >[V,,-Cos(@,t) + V,, - Sin(a,t)] 

k=0 
and by (carefully!) choosing sufficient irregularly-spaced time samples, 
the RHS coefficients may be determined. 
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Frequency-Domain Nonlinear Analysis Methods / 
Spectral Balance 


¢ Many of the problems with HB (especially in the multi-frequency case) 
are related to the FFT operation linking time- and frequency-domains. 
However, the FFT may be eliminated entirely by treating both linear and 
nonlinear parts of the system in the frequency domain 


¢ This is achieved by expressing signal variables at the interface by pre- 
defined truncated summations of complex phasors, and expressing the 
non-linear relationships between these in functional form; 

¢ Various forms of functional relationship have been described, including 
Generalised Power Series (Steer), Chebyshev expansions (Narhi), 
rational functions (Limiti)... ; 

¢ Practical application requires careful ‘bookkeeping’ to keep track of 
frequency components. The most complex operations are the nonlinear 
ones of multiplication/division (convolution/de-convolution); 
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Advantages of Harmonic Balance 

¢ Complex linear circuit environments may be easily and efficiently included 
in analysis; 

¢ Periodic solution can be obtained without needing a (perhaps) very-long 
transient period to elapse; 


* Solution can be very fast if non-linearity is moderate and/or a good initial 
guess is available. Recent improvements to the algorithm incorporating 
Krylov-subspace techniques have greatly extended the scale and 
complexity of problems which can be successfully tackled; 


Disadvantages of Harmonic Balance 
¢ Prior assumption of unique periodic solution may not always be justified; 


¢ In spite of great advances in numerical implementation, convergence 
failure (or extremely slow convergence) is still not uncommon at very large 
drive levels; 

e The method does not work well where the response of interest is not 
essentially steady-state but is instead a true transient. 
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Circuit Envelope Techniques 


¢ Several research groups have sought to combine the best features of 
Hamonic Balance and Time-Domain within a single simulator technology. 


¢ The Circuit Envelope simulator promoted by HP-EEsof, may be viewed as 
an extension to Harmonic Balance. Now each of the Fourier Series 
coefficients is considered as a complex-valued function of time which can 
represent an arbitrary modulating spectrum around each tone: 


+K : 
| u(t)= >U,(t)- eo 
| k=-K 


¢ A time-domain solution is performed consisting of a sequence of HB 
solutions at each time-step. However, the most critica! point is that the 
transient analysis time-step is determined by the modulation waveform, 
not by the high-frequency tone(s). Some modifications to standard HB are 
then required. For example, a linear capacitor now returns a k‘ harmonic 
current component given by: 


transient envelope 


dt 
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Diagrammatic Representation of Circuit Envelope 


Input Signal Circuit Envelope Output, = 


Amplitude & 
Phase ve g ane 


+ Amplitude ys. Time ¢ Spectral 
Time Phase vs. Time regrowth 
+ Frequency (S51) « Mixer 
vs. Time intermogs 
+ Harmonics vs. Time « Phase neise 
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Discrete Convolution Analysis 


In principle, convolution methods provide an attractive approach to 
incorporate arbitrary frequency domain blocks within a time-domain 
transient simulation. In practice, considerable care is needed in 
numerical implementations to achieve causality, stability and accuracy; 


Frequency-Domain: 


ad ()=¥(f)-V(A) 


Time-Domain: 


i(t) = Jay(t —t)-v(t)-dr 


Linear 1-Port where: h,(t)= [Y(f)-e/2*". df 
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-4GHZ -2GHz 0GHz +2GHz2z +4GHz 


Periodic System Function 
in Frequency-Domain has 
Exact Discrete Impulse 
Response Equivalence 


Impulse Weight 
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Two-Tone Intermodulation Test, PHEMT Amplifier 
f,= 10GHz Af=2.5MHz (Convolution Analysis) 
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Test with 8 Tones [f,= 10GHz, Af=2.5MHz] 


Load Voltage (V) 
Output Power (dBm) 


-80 
9.9 9.92 9.94 9,96 9.98 10 10.02 10.04 10.06 10.08 10.1 
Frequency (GHz) 
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Conclusions 


Modern high-frequency CAD provides a powerful set of tools 
which can provide key strategic advantages in meeting new 
design challenges; 

The down-side is that the range and sophistication of these tools 
is now placing considerable demands on the user; 

Nonlinear active device modelling remains a critical concern. 
Coping with emerging new technologies and process variations 
within existing technologies compound the modelling problem; 
Significant progress has been made in recent years in meeting 
the simulation challenges posed by complex RF systems 
carrying digitally modulated signals. 


University College Dublin 


4 , gy Tore inte e ba 


ai 
Oe eit eed te ioz J reuse “its en ‘i 
won Qoiieeny rl eegmityev ii op sonia 
alsiot aeortt te otal pits dprusy ert ‘ait ai obit ott. 
be ie air att, hd atigmenb | olde vehieacs srt inal WMO BF 
~) moones ipoltisd ganieiet oll aborn oaiveh av\tod penta 6 ae 
Zaitsev sae oie eh, re ati pata np Vhien tharos se a ig ietiqod fi ay | 
magico prt Hotionn en Grivocmny: sedpolonsts Gh tee cai | 
baifeem mlz ast miasenenten ‘rooad aur atorgac, preme e(P9 as 
ete vt 1x colanna9, ct terol polignti-1/4d'sh nena, 
=) eee 1s ANBIe poteidbar (On Tae en His 


x4 


. a 
2 4 * 
7 iP A 
L. * | i eee 
fal y ‘ . f P 
ere Seer zr ony “ane orn hart beeaia nc tpeme, Cunikplant "Greer NHN —etRenhmme ety tm smi aims mw S emmmetpeomene jatyalllnly ON ee ee tag 
ceiheyut), Bet le 2» er a) shiny id ae | 0 ice vaient 
ie ri ; ’ 
. ' 
- - ~~ mie a ~ _ ny gt ene comet ‘ 
4 
i ry 
1 ;¥ 
| a, 
' c ; x < 
A> “y i f 
_— at o tar ot - ° - os ee m en er 
5 
*4 
j y my 4 Ne 7 ‘ ‘ait & a i 
7 > U oe ’ , 2 } ie io ‘ 
ee, * i " ; ‘ee ‘ “4 ‘ Pre wh @ Pas) AY 
\ an FN ss ie pc ie t Oyo“ dean ths 1 0 +2, sh iy y 
ua Jeol (MRS TONGS Gt sf in Ps | ook 
rae 7 : 3 rr J = ia ete 
> 1 4 UG 
50 is ‘ hs 


Aare item ap iat 


ee | . 7 
oh ais 
i 


ia 3 rae 


= Aegan 


